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Abstract 

Pyrroloquinoline quinone (PQQ) is a redox cofactor utilized by a number of prokaryotic dehydrogenases. Not all prokaryotic 
organisms are capable of synthesizing PQQ, even though it plays important roles in the growth and development of many 
organisms, including humans. The existence of PQQ-dependent enzymes in eukaryotes has been suggested based on 
homology studies or the presence of PQQ-binding motifs, but there has been no evidence that such enzymes utilize PQQ as 
a redox cofactor. However, during our studies of hemoproteins, we fortuitously discovered a novel PQQ-dependent sugar 
oxidoreductase in a mushroom, the basidiomycete Coprinopsis cinerea. The enzyme protein has a signal peptide for 
extracellular secretion and a domain for adsorption on cellulose, in addition to the PQQ-dependent sugar dehydrogenase 
and cytochrome domains. Although this enzyme shows low amino acid sequence homology with known PQQ-dependent 
enzymes, it strongly binds PQQ and shows PQQ-dependent activity. BLAST search uncovered the existence of many genes 
encoding homologous proteins in bacteria, archaea, amoebozoa, and fungi, and phylogenetic analysis suggested that these 
quinoproteins may be members of a new family that is widely distributed not only in prokaryotes, but also in eukaryotes. 
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Introduction 

Pyrroloquinoline quinone (Figure 1, PQQ), a cofactor of 
prokaryotic dehydrogenases, was discovered as a redox cofactor 
of bacterial glucose dehydrogenase in 1964 [1] Its chemical 
structure was reported after its isolation from primary alcohol 
dehydrogenase of methylotrophic bacteria [2] . Some prokaryotic 
organisms are able to synthesize PQQ, whereas others require an 
exogenous source. The enzyme group to which the PQQ cofactor 
belongs is the quinoprotein group [3-5], and all PQQ-dependent 
quinoproteins reported to date have a characteristic propeller-fold 
superbarrel structure [3,5]. Bacterial PQQ-dependent dehydroge- 
nases which have an eight-bladed structure contain the charac- 
teristic amino acid sequence, which has been usually used to 
identify PQQ-dependent proteins. Although PQQ is physiologi- 



cally important in plants [6] as well as mammals (including 
humans) [7-10], there has been no evidence showing PQQ as a 
cofactor of any eukaryotic enzyme so far. 

A new class, named "Auxiliary Activities" (AA), was recently 
introduced in the database of carbohydrate-active enzymes 
(CAZy) to cover the full extent of the lignocellulose breakdown 
machinery [11]. This class covers a broad group of catalytic 
modules associated with plant cell wall degradation. The members 
of the AA family 8 consist of a cytochrome domain of spectral class 
b (protoheme IX); the first example of this family was found as a 
domain of cellobiose dehydrogenase (CDH; EC 1.1.99.18) isolated 
from the wood-rotting basidiomycete Phanerochaete chrysospor- 
ium, in which the cytochrome domain was linked to a flavin- 
containing AA family 3 flavin-containing sugar dehydrogenase 
domain [12]. We are interested in the AA family 8 hemoproteins 
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Figure 1. Chemical structure of PQQ. 
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because the heme domain is the first cytochrome ^-type heme 
bound through Met/His with a primarily (3-structure [13], and 
appears to be unrelated to any known protein other than the 
cytochrome domain of CDHs. 

The basidiomycete Coprinopsis cinerea (formerly known as 
Coprinus cinereus, Figure 2) produces many extracellular oxido- 
reductases, including peroxidases and laccases, as well as 
cellulolytic enzymes, and has many putative oxidoreductase genes. 
To identify novel proteins with this unique heme ^-containing 
cytochrome domain, the amino acid sequence of the cytochrome 
domain of CDH from P. chrysosporium was subjected to BLAST 
search using the C. cinerea Okayama 7 (# 1 30) genome database 
at the Broad Institute. We fortuitously discovered a novel PQQ; 
dependent protein that belongs to a category that is different from 
the known PQQ;quinoprotein family during a search for AA 
family 8 cytochrome domains. 

Materials and Methods 

Materials 

Coprinopsis cinerea strain 5338 was kindly provided by Dr. 
Yasuhiro Ito (National Food Research Institute, Ibaraki, Japan) 
[14]. Pyrroloquinoline quinone was purchased from Sigma- 




Figure 2. Fruiting body of Coprinopsis cinerea, 
doi:1 0.1 371/journal. pone.01 04851 .g002 



Aldrich (Tokyo, Japan). Bovine heart cytochrome c was purchased 
from Wako Pure Chemical Industries (Osaka, Japan). 

Cloning and transcript analysis 

C. cinerea was grown in Kremer and Wood medium [15] 
containing 2% Avicel (Merck, Whitehouse Station, NJ) for 12 d, 
and then total RNA was extracted using the RNeasy Plant Mini 
Kit (Qiagen, Hilden, Germany). First-strand cDNA was synthe- 
sized using Superscript III reverse transcriptase (Invitrogen, 
Carlsbad, CA) with an oligo(dT) primer (Takara Bio, Shiga, 
Japan). PCR using KOD plus version 2 DNA polymerase 
(Toyobo, Osaka, Japan) with the synthesized first-strand cDNA 
as a template was primed with oligonucleotides flanking the 
predicted translational start and stop codons (20 nucleotides 
upstream, 5'-TCGGGACCGACCACGAACG-3'; 63 nucleotides 
downstream, 5'-CGATTCTGTCTTGAAGCCCGACT-3'). The 
PCR products were subcloned into the pGEM-T Easy vector 
(Promega, Madison, WI) followed by sequencing with a 3130 
Genetic Analyzer (Applied Biosystems, Foster City, CA). The 
nucleotide sequence reported in this paper has been deposited in 
the DDBJ database under accession number AB901366. 

Sequence analysis 

The C cinerea genome database (http://www.broadinstitute. 
org/ annotation/ genome/ coprinus_cinereus/MultiHome.html) 
was searched with the amino acid sequence of CDH from P. 
chrysosporium (accession no. AAB92262) using the TBLASTN 
algorithm [16,17]. The search was carried out with standard 
settings and the BLOSUM 62 matrix. The amino acid sequences 
of candidate genes were scanned for the presence of signal peptides 
using SignalP version 4. 1 software [18] at the Center for Biological 
Sequence Analysis (http://www.cbs.dtu.dk/services/SignalP/) 
and searched against the protein database at the National Center 
for Biotechnology Information (http://blast.ncbi.nlm.nih.gov/ 
Blast. cgi). The homology structure model was constructed using 
Protein Homology/analogy Recognition Engine (Phyre) version 
2.0 (http://www.sbg.bio.ic.ac.uk/phyre2/) [19]. The CcSDH 
model was used for preparing a figure in which PQQ^ of sGDH 
from Acinetohacter calcoaceticus and heme in the cytochrome 
domain of P. chrysosporium CDH (PDBID ICRU and 1D7C, 
respectively) were superimposed to the model. The structure based 
alignment was created during the process of homology modeling 
using the Phyre2 server, and the corresponding figure was 
prepared using ESPript ver. 2.3 [20]. 

Production and purification of recombinant protein in 
P. pastoris 

The expression vector was constructed as described [21,22]. 
Two different oligonucleotide primers were designed based on the 
nucleotide sequence of the mature protein to allow for ligation into 
the EcoKl and Notl sites of the Pichia expression vector pPICZcx 
(Invitrogen, Carlsbad, CA): 5' primer (5'-TTGAATTC- 
CAAGGCTCTCCCACTCAGTG-3') and 3' primer (5'- 
TTGCGGCCGCTATGCAGGAACACACTGAGAGTACC- 
3'). PCR was carried out using the primer pair with KOD plus 
version 2 DNA polymerase (Toyobo, Osaka, Japan) and the first- 
strand cDNA synthesized above as a template. The PCR product 
was subcloned into the pGEM-T Easy Vector (Promega, Madison, 
WI) followed by sequencing with a 3130 Genetic Analyzer 
(Applied Biosystems, Foster City, CA). The target gene was 
digested with EcoKl and Notl and ligated into the corresponding 
restriction sites of the pPICZcxa vector. The nucleotide sequence of 
the inserted cDNA was confirmed by sequence analysis. Approx- 
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Figure 3. Domain organization of CcSDH, CDH from P. chrysosporium, CBMl -carrying CDH from T. heterothallica, and cellulose- 
binding cytochrome bsei from P. chrysosporium. Abbreviations: AA3, Auxiliary Activities (AA) family 3 enzymes defined as flavoproteins 
containing a flavin-adenine dinucleotide (FAD)-binding domain; AA8, AA family 8 enzymes with the cytochrome domain of spectral class b; CBMl, 
family 1 carbohydrate-binding module. 
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imately 10 |ig of the DNA construct in pPICZcx was linearized 
with Bpu 11021 (Takara, Japan) prior to transformation into P. 
pastoris. Electroporation, selection of the transformant, and 
production of recombinant protein were carried out using an 
EasySelect Pichia expression kit (version G; Invitrogen, Carlsbad, 
GA) according to the manufacturer's instruction. The PQQ^ 
domain of the recombinant protein was constructed by deleting 
the coding regions for the heme domain and the carbohydrate- 
binding module (GBM) domain of the recombinant protein by 
PGR using pPIGZcx containing the full-length gene as a template. 
The two oligonucleotide primers used to truncate the heme 
domain (221 amino acids from the N- terminus) were 5'- 



GGAGAAAAGAAGGTTGGTGTGTTGG-3' as the 5' primer 
and 5'-GGAAGAGAGGAAGGTTGTTTTGTGG-3' as the 3' 
primer; these primers included the (X-factor signal sequence of the 
pPIGZoc vector upstream of the recombinant protein gene and the 
gene coding for the amino acids after residue T222. For truncation 
of the GBM domain, a stop codon was introduced by using the 
following primers: 5'- GATGATGAAGGGGTAGTGGGGGGG- 
TATTGTTGAGG-3' (5' primer) and 5'- GGTGAAGAA- 
TAGGGGGGGAGTAGGGGTTGATGATG-3' (3' primer). 
The underlined codons represent mismatches that will introduce 
a stop codon and prevent the translation of the last 77 G -terminal 
amino acids, including the GBM domain. The mutations were 
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Figure 4. A, SDS-PAGE of purified (lane 1) and deglycosylated (lane 2) CcSDH. B, Effect of PQQ on cytochrome c-reducing activity of CcSDH 
with L-fucose as a substrate. For the activity measurement in B, 40 nM CcSDH was incubated with 20 mM L-fucose, 50 [iN\ cytochrome c and 1 |aM 
PQQ in HEPES buffer (pH 7.0), and the reduction of cytochrome c was monitored based upon the increase of absorption at 550 nm. C and D, 
Isothermal titration calorimetry analysis of the PQQ domain of CcSDH and PQQ with the raw data (C) and plots of the integrated peaks (D). The PQQ 
domain of CcSDH (6.73 laM) was titrated with an initial injection of 2 [lL of PQQ (46.3 iiM) followed by 69 consecutive injections of 4 [lL of PQQ in the 
presence of 1 mM CaCb and 20 mM sodium acetate (pH 6.0). The black line indicates the best fit to a single-site model. 
doi:1 0.1 371/Journal.pone.01 04851 .g004 
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D-arabinose 
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L-gulose 
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0.63 


D-lyxose 


12.9 (±0.3) 


66.8 (±2.6) 


0.19 
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confirmed by DNA sequencing, and the recombinant plasmid was 
transformed into P. pastoris for subsequent protein expression 
according to the same protocol used for the full-length recombi- 
nant protein. 

The culture was centrifuged (30 min, 1,500 xg), and the crude 
proteins (including the recombinant protein) were precipitated 
with ammonium sulfate (70% saturation) from the cell-free culture. 
The precipitates were dissolved in 20 mM Tris-HCl buffer 
containing 1 M ammonium sulfate (pH 8.0). The crude proteins 
were fractioned on a Toyopearl Phenyl-650S column equilibrated 
with 20 mM Tris-HCl buffer containing 1 M ammonium sulfate 
(pH 8.0). The recombinant protein was eluted using a linear 
reverse gradient to 20 mM Tris-HCl buffer (pH 8.0). The fraction 
containing the recombinant protein was collected and equilibrated 
against 20 mM Tris-HCl buffer (pH 8.5). The solution was loaded 
on a Toyopearl DEAE-650S column equilibrated with 20 mM 
Tris-HCl (pH 8.5). The enzyme was eluted from the column using 
a linear gradient in the same buffer (0 to 500 mM NaCl). The 
purity was confirmed by SDS-PAGE analysis on 12% polyacryl- 
amide gels, and the N-terminal amino acid sequence was 
determined with a protein sequencer (model 491 cLC; Applied 
Biosystems, Foster City, CA) as described. The recombinant 
protein concentration in the solution was estimated from the 
absorbance at 420 nm (8420= 130 mM Vm ^). Deglycosylation 
of the recombinant protein was performed using endo-P-A^- 
acetylglucosaminidase H (endo-H) as described [21,22]. 

Enzynne activity 

The enzyme activity was assayed using cytochrome c as an 
electron acceptor [23,24]. The assay was performed by photo- 
metric monitoring of the time-dependent reduction of 50 |iM 
cytochrome c at 550 nm (8550= 17.5 mM ^ cm ^) at 30°C. One 
activity unit (U) corresponds to the amount of enzyme that can 
convert 1 |imol of substrate per minute. The substrate specificity 
of the recombinant protein (40 nM) was assessed using 1 mM D-, 
L-allose, D-, L-galactose, D-, L-glucose, D-, L-gulose, D-, L- 
mannose, D-, L-talose, D-, L-arabinose, D-, L-lyxose, D-, L- 
xylose, D-, L-fructose, D-, L-tagatose, L-sorbose, cellobiose, 
maltose, sophorose, trehalose, D-, L-fucose, L-rhamnnose, N- 
acetyl-D-glucosamine, D-glucosamine, and D-glucosone. To 
obtain the kinetic parameters of oxidation for L-gulose, D- 
arabinose, D-lyxose, L-fucose, and D-glucosone, the assay was 
performed with various concentrations of the saccharides toward 
50 nM recombinant protein in 50 mM Tris-HCl buffer (pH 8.5) 
containing 1 |iM PQQ^ and 1 mM CaCl2. To determine the 
Michaelis constant {K^) and catalytic rate (^cat)? the experimental 
data were fitted to the Michaelis-Menten equation. 



Isothernnal titration calorimetry (ITC) 

ITC experiments were carried out at 25°C with a MicroCal VP- 
ITC (GE Healthcare, Northampton, MA, USA). A sample of the 
purified PQQ^ domain of the recombinant protein was prepared by 
dialysis against ITC buffer (20 mM sodium acetate, pH 6.0, 
1 mM CaCl2). The concentration of the PQQ^ domain of the 
recombinant protein was determined from the UV absorbance at 
280 nm using the calculated molar extinction coefficient 
7.241 X 10* M"Vm"\ PQa was dissolved in the same ITC 
buffer, and the concentration was measured based upon the 
absorbance at 257 nm using the molar extinction coefficient [25] 
2.15xlO^M" cm . Solutions were de-gassed for 1 min using a 
vacuum degasser immediately prior to experimentation. The 
titration sequence consisted of an initial injection of 2 |iL of PQQ^ 
(46.3 |iM) followed by 69 additional injections of 4 |iL of PQQ^at 
220 s intervals into a calorimeter cell containing the PQQ^ domain 
of the recombinant protein solution (6.73 |iM) stirred at 300 rpm. 
Baseline-subtracted data were fitted according to Marquardt's 
single-site model using Origin ITC Analysis software (Microcal 
Software, Northampton, MA) to obtain stoichiometric and 
thermodynamic parameters. 

Construction of a phylogenetic tree 

The National Center for Biotechnology Information protein 
database was searched for amino acid sequences corresponding to 
the catalytic domain of the recombinant protein using the PSI- 
BLAST algorithm. All searches were performed with standard 
settings and the BLOSUM 62 matrix. The sequences obtained 
and the sequences of known quinoproteins (PQQ;dependent 
methanol dehydrogenases, ethanol dehydrogenases, soluble glu- 
cose dehydrogenases, quinohemoprotein alcohol dehydrogenases, 
polyvinyl alcohol dehydrogenase, membrane-bound glucose 
dehydrogenase, membrane-bound alcohol dehydrogenase, and 
aldose sugar dehydrogenases) were subjected to alignment analysis 
using MAFFT (ver. 6.85) [26,27]. Non-conserved regions were 
trimmed using SeaView (ver. 4.4.2) software. The phylogenetic 
tree was generated from the trimmed sequences using ClustalX 
(ver. 2.1) with the neighbor-joining method (Bootstrap value, 
1000). 

Results 

As a result of the BLAST search, several genes, including genes 
homologous to the cytochrome domain of CDH, were identified in 
the C. cinerea genome, and we cloned one of them (Chromosome 
6:1274274-1277156). RT-PCR using total RNA of C. cinerea was 
performed to obtain the cDNA of the whole protein. The cDNA 
consisted of 2883 bp, which included an open reading frame for 
726 amino acids. Based on the analysis of the predicted amino acid 
sequence using the SignalP 4.1 server [18], the N-terminal 18 
amino acid sequence was identified as a signal peptide, suggesting 
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Figure 5. Phylogenetic tree of prokaryotic and eukaryotic quinoproteins. The phylogenetic tree was generated from the amino acid 
sequences of CcSDH, proteins homologous to CcSDH, and known quinoproteins using ClustalX (ver. 2.1) with the neighbor-joining method. 
doi:1 0.1 371 /journal. pone.01 04851 .g005 



that the protein is secreted extracellularly. Homology searches of 
the amino acid sequence in the NCBI protein database revealed 
that the protein was composed of an N-terminal AA family 8 
cytochrome domain and a C -terminal cellulose-binding domain 
belonging to the CBM family 1 as shown in Figure 3. In addition 
to these domains, the protein contains an unknown domain, which 
has a low similarity to glucose/sorbosone dehydrogenases, 
different than the AA family 3 flavin domain seen in CDHs. 

The recombinant protein was named C. cinerea sugar 
dehydrogenase (CcSDH) because of its oxidative activity towards 
various sugars. This enzyme was heterologously expressed in the 



methylotrophic yeast Pichia pastoris and purified from the culture 
solution using a two-step column chromatography. The purified 
enzyme yielded a single band at a molecular mass of 98 ±8 kDa on 
SDS-PAGE (Figure 4A). The molecular mass measured was larger 
than that estimated from the amino acid sequence (77 kDa), 
suggesting that CcSDH is A^- and 0-glycosylated, as predicted by 
the NetNGlyc and NetOGlyc servers. This conclusion was 
confirmed by a decrease of the apparent molecular weight 
following treatment with endo-glycosidase H as shown in 
Figure 4A. In anticipation that the unknown domain would have 
redox activity based on the electron transfer ability of the 
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Figure 6. 3D structure model of CcSDH showing its three domains the N-terminal AA family 8 cytochrome domain (orange), the 
PQQ-dependent dehydrogenase domain (magenta), and the C-terminal cellulose-binding domain belonging to CBM family 1 

(green), PQQ of sGDH from Acinetobacter calcoaceticus and heme in the cytochrome domain of P. chiysosporium CDH (PDBID 1CRU and 1D7C, 
respectively) were superimposed to the model. 
doi:1 0.1 371 /journal. pone.01 04851 .g006 



cytochrome domain, we tested NAD(P), FAD, and PQQ^ as 
cofactors for two-electron oxidation reaction of various sugars, 
alcohols, aldehydes, and their derivatives. The assay was carried 
out by monitoring the rate of cytochrome c reduction (similar to 
the CDH assay [23,24]). Interestingly, when various sugars 
(including D-glucosone, L-fucose, and some rare sugars) were 
used as substrates, the reduction of cytochrome c was observed 
only in the presence of PQQ. CcSDH had catalytic activity toward 
L-galactose, L-gulose, D-talose, D-arabinose, D-lyxose, L-fucose, 
and D-glucosone in the presence of PQQ and CaGl2. CcSDH 
shows little or no activity toward abundant sugars such as D- 
glucose, D-fructose, and cellobiose. The steady-state kinetic 
parameters for L-gulose, D-arabinose, D-lyxose, L-fucose, and 
D-glucosone were measured using the cytochrome c assay under 
optimal condition as summarized in Table 1. Because D- 
glucosone has the reductive activity of cytochrome c, we 
determined the initial velocity in an enzymatic reaction after 
subtracting the rate of cytochrome c reduction by D-glucosone 
itself from an apparent rate of cytochrome c reduction. The 
catalytic efficiency measured was in the following order: D- 
glucosone>L-fucose>D-arabinose>L-gulose>D-lyxose. D-gluco- 
sone exhibited the highest catalytic efficiency, k^at^K^, calculated 
to be 9.37 x 10^ s~^M~\ The values for these sugars appear to 
be large. However, considering the fact that the PQQ-dependent 
sugar dehydrogenases, such as soluble glucose dehydrogenase from 
Acinetobacter calcoaceticus (sGDH) have a large value 
(25 mM) for glucose, a physiological substrate [28], the value 
for D-glucosone (7.9 mM) with CcSDH is reasonable for a 
quinoprotein. Figure 4B shows the time-course of the absorbance 
at 550 nm based upon the reduced cytochrome c for L-fucose as a 
typical example. When CcSDH was added to the assay solution 
without PQQ no catalytic activity was observed. Upon addition of 
PQQ the absorbance immediately increased. To determine the 
binding constant of the protein to PQQ the PQQ domain alone 
was expressed in the Pichia expression system and the binding 
constant was estimated by ITC as shown in Figure 4C. The 
dissociation constant {K^ was 1.11 nM, and the molar ratio of 
PQQ to the domain was 1:1 (Figure 4D). 



BLAST search revealed many genes that encode enzymes 
similar to CcSDH in bacteria, archaea, amoebozoa, and fungi. 
Therefore, this type of quinoprotein appears to be widely 
distributed not only in prokaryotes, but also in eukaryotes. A 
phylogenetic tree generated from the amino acid sequences 
showed that these quinoproteins are distinct from other classes 
of known quinoproteins (Figure 5), suggesting that they represent a 
new family. 

Discussion 

The highest amino acid homology to the unknown domain of 
CcSDH was obtained with Deinococcus radiodurans L-sorbosone 
dehydrogenase (26%). Glucose dehydrogenase and sorbosone 
dehydrogenase, which have the highest degree of similarity to this 
protein, are both PQQ-dependent dehydrogenases. The catalytic 
activity for various sugars was observed only in the presence of 
PQQ. The strong binding activity of the protein to PQQ was 
shown by ITC. Based on these findings, we concluded that 
CcSDH is a PQQ-dependent enzyme, even though it shows low 
amino acid sequence homology with known PQQ-dependent 
enzymes. 

The three-dimensional structures of PQQ-dependent quino- 
proteins are generally divided into two types: an eight-bladed beta- 
propeller structure with each blade consisting of a four-stranded 
anti-parallel beta-sheet, and a six-bladed beta-propeller structure 
[3,29]. These two types of enzymes have no amino acid sequence 
homology with each other. The three-dimensional structure of 
CcSDH was modeled using the Phyre2 protein fold recognition 
server and the predicted structure was compared with those of 
other quinoproteins and structurally similar proteins. CcSDH was 
predicted to have a six-bladed instead of an eight-bladed 
quinoprotein structure (Figure 6). As shown in Table 2, CcSDH 
shows some structural homology with known prokaryote PQQ- 
dependent sugar dehydrogenases, sGDH [30-32] and soluble 
aldose sugar dehydrogenases (Asd) from Pyrobaculum aerophilum 
[33] and Escherichia coli [34], although the amino acid sequence 
homology is low (—20%). Among prokaryotic sugar dehydroge- 
nases, CcSDH shows some structural similarity with human 
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VGALHSGDGSQRLFILEKEGYVKILTPEGEIFKEPYLDIHKLVQSG . . . .IK 

HALLWGPDN. . Q I WL T E R AT G K I L R VN P E S G S V K T V F. . . .QVPEIVN 

WSIAPLGGG. . RYLVTERPGRLVLI SPSGKKLVASF D. . . . VA 

WALAFLPDN. . H GML I T L R G G E L R H W Q AG K G L S AP L S G V P D. . . . VW 

WGLAPLPGG. .DLLVSSRDEATITRVDAKTGRKTEL GEVPGVS 



CcSDH 

2WG3 

ICRU 

3A9G 

2G8S 

3DAS 



310 
278 
72 
86 
59 
76 



. IQDTQINHGIDVH 
. GGDERGLLSLAFH 
DADGQNGLLGFAFH 
. NVGEAGLLGLALH 
. AHGQGGLLDVVLA 
. PSGEGGLLGIALS 



S....G...R RIIA S SG DIAWSW DYDPATM 

NYKK.N. . .G KLYV S Y T T N Q E R W A I G P H D H I LRVVEY TVSRKNP 

DFKN.N. . .P YIYISGTFKNPKSKELP NQ TIIRRY TYNKSTD 

EFPKKS. . .W VYLY A S YF AE GGH I RNRV I RG RLDGSTF 

D. . . .FAQSR RIWL SYSEVGDDGK AG TAVGYG RLSDDLS 

D. . . .Y. . .ASDHMVYAYF T SA S D N R I VRM L Y D E K K P S G E Q 



CcSDH 347 T . . . . ATNRRTL . VTG . . . MNNF Y . 

2WG3 332 HQVDLRTARVFL . E VA . . . ELHRK . 

ICRU 127 T . . . . LEKPVDLLAGL . . . PSSKD . 

3A96 135 K . . . . LKEVKTL . IDG . . . IPGAYI 

268S 108 K . . . . VTDFRTV . FRQMPKLSTGN . 

3DAS 124 L . . . . GAPDTVF . RGI . . . PKGVI . 



FTRTVHISRKYPNLFALNVGSDG.NIDVPTRQQ NSGRAQIRVGDYDQ. 

LGGQLLFGPDGFLYIILGDGMITLDDMEEMDGL SDFTGSVLRLDVDTD 

QSGRLVIGPDQKIYYTIGDQGRN.QLAYLFLPNQAQHTPTQQELNGKDYHTYMGKVLRLNLDG. 

NGGRIRFGPDGMLYITTG . . . D A . A D P R L AQ D L S S L AG K I L R V D E E G . 

FGGRLVFDGKGYLFIALGENNQR.PTAQDLDKL QGKLVRLTDQGEIP. 

NGGRI AFGPDKMLYAGTGE SGD . . . .TGLSQDR K S L G GK I L RM T P D G . 



CcSDH 

2WG3 

ICRU 

3A96 

2G8S 

3 DAS 



410 LPQNGVP . FVSQYG . . . R VL G yHl RN D VG I T E D R . AG N I H S I E N S L D N A Y RM VN G Q R R D I H T N N P A E K V Y N L G D P S N RR A I F GG 



401 MCNVPYSIPRSNP . HFNSTNQPPEVFAH 

207 SIPKDNPSFNGVVS . . . HI YTL 

196 RPPADNP . FP . . NS . 

174 DDNPFIK . ESGVRA . 



184 EPAP . . G . NPFPGS . 



. PIWSY 
. EIWSY 
. P VYS Y 



LHDPGRCAVDR.HPTDI NINLT ILCSDSS 

HRNPQGLAFTP.NGKLLQSEQGPNS DDEINLIVKGGNYG 

HRNPQGIDWHRASGVMVATEHGPVGHDEVNIILK GGNYG 

IRNPQGMAMNPWSNALWLNEHGPRG GDEINIPQKGKNYG 

HRMVQGLAWDD .KQRLFASEFGQDTWDELNAIKP GDNYG 

AA 



CcSDH 489 YPDCYTVWEPSDFTDSPKQPGDWFTQDNSGQYTDAWCNANAVKPTLLL PPHTAPLDM 

2W63 463 ARILQI I YESEPSLL EFKPPLVGG 

ICRU 265 WPNVAGYKDDSGYAYAN YSAAANKS I KDLAQNGVKVAAGVPVTKE SEWTGKNFVPPLKTLYTVQDTYNYNDPTCGEMT Y ICWPTVAPS SA 

3A9G 252 WP L AT G K AG RGE F V D P V I D T G SETWAPSGA 

2G8S 232 WPLATWGINYSGFKIPEAKG E I V AG T E Q P V F Y W KDSPAVSGM 

3DAS 239 WP E AE G K GG G S G F H D P V AQ W STDEASPSG 



CcSDH 546 KFGLGN . DT N L Y VAL H G S WNRQP P QG YK VV WP GQ Y S AS G E W S P T AP L AQ S RT AW S D L LTN RNE N Q C S G F GN AN C F RP VG L VW S ADGQ N L 

2WG3 495 FVYRGC.QSE RLYGSYVFGDRN.GNFLTLQQQWQEKPLCLGTSGSCRGYFSGHILGFGEDELGE.V 

ICRU 355 YVYKGGKKAIT. . G WE N T L L VP S LPR G V I F R I K L D P T Y S T T Y D D A V P M F K S N N R Y RD V I A S P D G N . V 

3A9G 282 SFVHGD.MFPG. . . L R G WL L I AC LRG S M L A A VN F G D N M E V RKISTFFKNVFGRLRDVVIDDDGG.I 

2G8S 274 AFYNSD.KFP. . . Q W Q Q KL F I G ALKD K D V I VM S V N G D K VT E D G R I L T D R G QR I RD VR T G P D G Y L Y 

3DAS 255 lAYAEG.SVWMAGLRGERLWRIPLKGT A A A AD P Q A F L E G E Y GRL RT V A P AG G D K . L 

A 

CcSDH 635 YVS SDT S GEVFIIKR 

2WG3 563 YILSSSKSMTQTH N GKLYKIVD 

ICQl 419 LYV LTD T AGNVQ . . . 

3A9H 343 LIS TSNRDGIGSLRAGD DKILKIV. 

268S 335 VLT DES S GELLKVS. 

3DAS 322 WLV TSN T D G^G D AK GG D D R I L E L E V 



Figure 7. Alignment of the amino acid sequences of the catalytic domain of CcSDH and known structure of six-bladed 
quinoproteins. Perfect matches are enclosed in boxes with a black background. Boxes with a magenta background indicate the amino acid 
residues interacting with PQQ via direct hydrogen bonds in known structures. The filled arrowhead indicates the position of the proposed catalytic 
histidine in bacterial quinoproteins. Conserved residues involved in PQQ-binding in known quinoproteins are indicated by open arrowheads if also 
conserved in CcSDH [32-34]. 
doi:1 0.1 371 /journal.pone.01 04851 .g007 



hedgehog-interacting protein (HHIP) [35]. As seen in the 
alignment of the amino acid sequence of CcSDH with those of 
known six-bladed quinoproteins (Figure 7), a putative catalytic 
histidine is conserved in all proteins, whereas several amino acids 
contributing to PQQ^ binding in the prokaryote sugar dehydro- 
genase are missing in CcSDH and HHIP. However, as demon- 
strated above, the activity of CcSDH is PQQ;dependent. Thus, 



the amino acids involved in the adsorption of PQQ^ should be 
different in this enzyme family, indicating the diversity of PQQ; 
binding motifs in eukaryotic enzymes. 

The PQQ;dependent sugar dehydrogenase that we discovered 
here has interesting enzymatic properties associated with charac- 
teristic cytochrome (AA family 8) and cellulose-binding domains 
(CBM family 1) as well as the enzymes (proteins) related to plant 



Table 2. Summary of sequence homology modeling results using the Phyre2 server. 





Enzyme 


Organism 


PDBID 


Identity 


Confidence 


E-value^ 


soluble glucose dehydrogenase (sGDH) 


Acinetobacter calcoaceticus 


ICRU 


16 


100.0 


1 .90E-34 


aldose sugar dehydrogenase 


Pyrobaculum aerophilum 


3A9G 


22 


100.0 


1.10E-28 


soluble aldose sugar dehydrogenase (Asd) 


Escherichia coli 


2G8S 


15 


100.0 


1 .30E-26 


human hedge hog interacting protein (HHIP) 


Homo sapiens 


2WG3 


16 


100.0 


1 .60E-27 


putative glucose dehydrogenase 


Thermus thermophilus 


2ISM 


16 


100.0 


9.00E-27 


aldose sugar dehydrogenase (Adh) 


Streptomyces coelicolor 


3DAS 


17 


100.0 


9.40E-24 



^E-value generated using the Phyre2 profile-profile alignment algorithm. 
doi:1 0.1 371/journal.pone.01 04851 .t002 
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cell wall degradation, and shows high activity towards rare sugars. 
Therefore, the discovery of the PQQjdependent domain in 
CcSDH could form the basis for a new AA family in the CAZy 
database. Detailed enzymatic studies will be necessary to evaluate 
why eukaryotes produce and secrete such enzymes, and how these 
enzymes acquire PQQ, because only a few prokaryotic organisms 
has been reported to be capable of synthesizing PQQ; Interest- 
ingly, genes encoding enzymes similar to CcSDH were found in 
various bacteria, whereas only a few archaea contained corre- 
sponding genes. Furthermore, all bacterial and archaeal CcSDH- 
like enzymes lack the cytochrome domain present in some of the 
fungal enzymes. This observation raises the question as to how 
eukaryotes such as fungi acquired genes encoding PQjQjdepen- 
dent enzymes. 

References 

1. Hauge JG (1964) Glucose Dehydrogenase of Bacterium Anitratum: An enzyme 
with a novel prosthetic group. J Biol Chem 239: 3630-3639. 

2. Salisbury SA, Forrest HS, Cruse WB, Kennard O (1979) A novel coenzyme 
from bacterial primary alcohol dehydrogenases. Nature 280: 843-844. 

3. Anthony G, Ghosh M (1998) The structure and function of the PQQjContaining 
quinoprotein dehydrogenases. Prog Biophys Mol Biol 69: 121. 

4. Anthony G (2001) Pyrroloquinoline quinone (PQQ) and quinoprotein enzymes. 
Antioxid Redox Signal 3: 757-774. 

5. Matsushita K, Toyama H, Yamada M, Adachi O (2002) Quinoproteins: 
structure, function, and bio technological applications. Appl Microbiol Biotech- 
nol 58: 13-22. 

6. Ghoi O, Kim J, Kim JG, Jeong Y, Moon JS, et al. (2008) Pyrroloquinoline 
quinone is a plant growth promotion factor produced by Pseudomonas 
fluorescens B16. Plant Physiol 146: 657-668. 

7. KiUgore J, Smidt C, Duich L, Romero-Ghapman N, Tinker D, et al. (1989) 
Nutritional importance of pyrroloquinoline quinone. Science 245: 850-852. 

8. Stites TE, MitcheU AE, Rucker RB (2000) Physiological importance of 
quinoenzymes and the O-quinone family of cofactors. J Nutr 130: 719-727. 

9. Misra HS, Rajpurohit YS, Khairnar NP (201 2) Pyrroloquinoline-quinone and its 
versatile roles in biological processes. J Biosci 37: 313-325. 

10. Harris GB, Ghowanadisai W, Mishchuk DO, Satre MA, Slupsky GM, et al. 
(2013) Dietary pyrroloquinoline quinone (PQQ) alters indicators of inflamma- 
tion and mitochondrial-related metabolism in human subjects. J Nutr Biochem 
24: 2076-2084. 

1 1 . Levasseur A, Drula E, Lombard V, Coutinho PM, Henrissat B (20 1 3) Expansion 
of the enzymatic repertoire of the CAZy database to integrate auxiliary redox 
enzymes. Biotechnol Biofuels 6: 41. 

12. Ayers AR, Ayers SB, Eriksson KE (1978) Cellobiose oxidase, purification and 
partial characterization of a hemoprotein from Sporotrichum pulverulentum. 
Eur J Biochem 90: 171 181. 

13. Hallberg BM, Bergfors T, Backbro K, Pettersson G, Henriksson G, et al. (2000) 
A new scaffold for binding haem in the cytochrome domain of the extracellular 
flavocytochrome cellobiose dehydrogenase. Structure 8: 79 88. 

14. Yanagi SO, Kawasumi T, Takebe I, Takemaru T (1988) Genetic analyses of 
Coprinus cinereus strains derived through intraspecific protoplast fusion. 
Agricultural and Biological Chemistry 52: 281-284. 

15. Kremer SM, Wood PM (1992) Evidence that cellobiose oxidase from 
Phanerochaete chrysosporium is primarily an Fe(III) reductase. Kinetic 
comparison with neutrophil NADPH oxidase and yeast flavocytochrome ^2- 
Eur J Biochem 205: 133 138. 

16. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local 
alignment search tool. J Mol Biol 215: 403 410. 

17. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped 
BLAST and PSLBLAST: a new generation of protein database search 
programs. Nucleic Acids Res 25: 3389-3402. 

18. Petersen TN, Brunak S, von Heijne G, Nielsen H (2011) SignalP 4.0: 
discriminating signal peptides from transmembrane regions. Nat Methods 8: 
785-786. 

19. Kelley LA, Sternberg MJ (2009) Protein structure prediction on the Web: a case 
study using the Phyre server. Nat Protoc 4: 363-371. 

20. Gouet P, CourceUe E, Stuart DI, Metoz F (1999) ESPript: analysis of multiple 
sequence alignments in PostScript. Bioinformatics 15: 305-308. 



It is worth noting that this discovery would not have been 
possible using the usual approach of searching for PQQ; 
dependent enzymes. Our results indicate the existence of a 
previously unknown enzyme family. That PQQ is a beneficial 
vitamin for humans has been reported [36]. In view of our 
findings, novel PQQjdependent enzyme(s) may exist in humans. 

Author Contributions 

Conceived and designed the experiments: MS HO MY KI NN. Performed 
the experiments: HM KU KT NS TL Analyzed the data: HM KT NS TI 
MY KI NN. Contributed reagents /materials /analysis tools: HM KU KT 
NS TI MY KI NN. Contributed to the writing of the manuscript: HM KT 
MY KI NN. 



21. Yoshida M, Igarashi K, Wada M, Kaneko S, Suzuki N, et al. (2005) 
Characterization of carbohydrate-binding cytochrome 6562 from the white-rot 
fungus Phanerochaete chrysosporium. Appl Environ Microbiol 71: 4548-4555. 

22. Yoshida M, Ohira T, Igarashi K, Nagasawa H, Aida K, et al. (2001) Production 
and characterization of recombinant Phanerochaete chrysosporium cellobiose 
dehydrogenase in the methylotrophic yeast Pichia pastoris. Biosci Biotechnol 
Biochem 65: 2050 2057. 

23. Samejima M, Eriksson KE (1992) A comparison of the catalytic properties of 
cellobiose: quinone oxidoreductase and cellobiose oxidase from Phanerochaete 
chrysosporium. Eur J Biochem 207: 103-107. 

24. Igarashi K, Verhagen MF, Samejima M, Schulein M, Eriksson KE, et al. (1999) 
Cellobiose dehydrogenase from the fungi Phanerochaete chrysosporium and 
Humicola insolens. A flavohemoprotein from Humicola insolens contains 6- 
hydroxy-FAD as the dominant active cofactor. J Biol Chem 274: 3338 3344. 

25. Jongejan JA (1989) Chemistry of PQQ. Delft, The Netherlands: Delft University 
of Technology. 

26. Katoh K, Misawa K, Kuma K, Miyata T (2002) MAFFT: a novel method for 
rapid multiple sequence alignment based on fast Fourier transform. Nucleic 
Acids Res 30: 3059-3066. 

27. Katoh K, Standley DM (2013) MAFFT multiple sequence alignment software 
version 7: improvements in performance and usability. Mol Biol Evol 30: 772- 
780. 

28. Igarashi S, Hirokawa T, Sode K (2004) Engineering PQQ glucose dehydroge- 
nase with improved substrate specificity. Site-directed mutagenesis studies on the 
active center of PQQ glucose dehydrogenase. Biomol Eng 21: 81 89. 

29. Toyama H, Mathews FS, Adachi O, Matsushita K (2004) Quinohemoprotein 
alcohol dehydrogenases: structure, function, and physiology. Arch Biochem 
Biophys 428: 10-21. 

30. Oubrie A, Rozeboom HJ, Kalk KH, Olsthoorn AJ, Duine JA, et al. (1999) 
Structure and mechanism of soluble quinoprotein glucose dehydrogenase. 
EMBOJ 18: 5187-5194. 

31. Oubrie A, Rozeboom HJ, Kalk KH, Duine JA, Dijkstra BW (1999) The 1.7 A 
crystal structure of the apo form of the soluble quinoprotein glucose 
dehydrogenase from Acinetobacter calcoaceticus reveals a novel internal 
conserved sequence repeat. J Mol Biol 289: 319-333. 

32. Oubrie A, Rozeboom HJ, Dijkstra BW (1999) Active-site structure of the soluble 
quinoprotein glucose dehydrogenase complexed with methylhydrazine: a 
covalent cofactor-inhibitor complex. Proc Natl Acad Sci U S A 96: 11787 
11791. 

33. Sakuraba H, Yokono K, Yoneda K, Watanabe A, Asada Y, et al. (2010) 
Catalytic properties and crystal structure of quinoprotein aldose sugar 
dehydrogenase from hyperthermophilic archaeon Pyrohaculum aerophilum. 
Arch Biochem Biophys 502: 81 88. 

34. Southall SM, Doel JJ, Richardson DJ, Oubrie A (2006) Soluble aldose sugar 
dehydrogenase from Escherichia coli: a highly exposed active site conferring 
broad substrate specificity. J Biol Chem 281: 30650-30659. 

35. Bishop B, Aricescu AR, Harlos K, O'Callaghan CA, Jones EY, et al. (2009) 
Structural insights into hedgehog ligand sequestration by the human hedgehog- 
interacting protein HHIP. Nat Struct Mol Biol 16: 698-703. 

36. Kasahara T, Kato T (2003) Nutritional biochemistry: A new redox-cofactor 
vitamin for mammals. Nature 422: 832. 



PLCS ONE I www.plosone.org 



8 



August 2014 | Volume 9 | Issue 8 | el04851 



